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Using an anisotropic hybridization matrix element, V(/c), reflecting cubic symme- 
try at the "rare-earth" sites of the infinite U Anderson lattice, we calculate the bare 
quasiparticle scattering amplitude to order 1/N using slave Bosons and find signifi- 
cant differences with previous results under spherical symmetry. Due to exchange of 
density fluctuations, we find an instability, toward E 9 and Ti s (even-parity) pairing. 
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Since the advent of heavy Fermion superconductivity as discovered in CeCu2Si2[l], physi- 
cists have puzzled over the cause and nature of the pairing in these materials. Unique for 
several reasons, heavy Fermion compounds are comprised of intermetallics and rare-earth or 
actinide atoms (such as U or Ce) characterized by: partially filled 4f or 5f shells; a large local 
Coulomb repulsion; and a rich multiplet structure in the presence of crystal electric fields. 
Hybridization between conduction electrons and these cerium or uranium electrons gives 
rise to extremely heavy quasiparticles which participate in the pairing [1]. These properties, 
coupled with the fact that the quasiparticle bandwidth (set by the lattice Kondo tempera- 
ture T ) is much less than the Debye temperature, make heavy Fermions strong candidates 
for non-phononic, anisotropic superconductors [2], [3]. Our work focuses on calculating the 
two-particle scattering amplitude for CeCu 2 Si 2 , and we see that the inclusion of crystal elec- 
tric fields of cubic symmetry at the Ce sites strongly modifies the anisotropy and strength 
of the quasiparticle interactions, leading to a superconducting instability of either E g or T ig 
symmetry. 

Most agree that a lattice of local electrons on "rare-earth" sites hybridized with con- 
duction electrons captures the essential physics of heavy Fermions in the normal state. 
Anderson suggested [4] that the inverse degeneracy of the "rare-earth" multiplet would be 
a useful expansion parameter to quantify the strength of quasiparticle interactions in the 
lattice. In the limit of infinite Coulomb repulsion, Coleman [5] made a low temperature 
(T<T G ), strong-coupling calculation tractable by introducing slave bosons, b(-Rj), (as first 
used by Barnes [6]) to represent a hole or empty site on the lattice. A physical electron on 
a rare-earth site is then created by the product of a fermionic operator and a slave boson: 
v ' r a (Ri)b(Ri) , where T and a are crystalline quantum numbers at site i. There is also a 
Lagrange multiplier to constrain the number of Fermions and Bosons at each site to unity. 

Lavagna, Millis, and P. Lee [7] looked for a superconducting instability in the Anderson 
lattice using this formalism. Within the so-called SU(N) version of the model, where one 
assumes N-fold degenerate conduction states, and with a spherically symmetric hybridization 
matrix element V, they found a d-wave pairing instability at order 1/N. Zhang and T. K. 

2 



Lee [8] pointed out the need to include spin-orbit coupling and an anisotropic hybridization 
in order to model heavy Fermion systems realistically. They showed that even-parity pairing 
channels had large local, or "hard-core", repulsions which prevented pairing. 

Our work reveals that the local repulsive interactions are substantially changed due to 
the cubic crystal fields, and this change is in large part responsible for the pairing instability. 
The result thus can be seen as a natural evolution of previous calculations [7,8], wherein the 
crystal field split multiplet structure is, in some sense, the last major piece of local physics 
to be added to the model. 

CeCu2Si2 has a superconducting transition temperature of about 0.6 K and a Kondo tem- 
perature (effective quasiparticle bandwidth) of approximately 10 K (T o ^10K) [1]. Inelastic 
neutron scattering measurements by Horn et al.[9] clearly show the crystal field splitting of 
the low-lying Ce if 1 spin-orbit multiplet (J— 5/2). Although CeCu 2 Si 2 has a tetragonal unit 
cell, fits to specific heat [10], static susceptibility [11], and de-Haas van-Alphan data for the 
effective mass[12] require that the Ce ions actually occupy sites of pseudo-cubic symmetry. 
Thus, our work focuses on plane wave conduction states hybridizing with localized crystal 
fields states of cubic symmetry. 

In the presence of crystal electric fields of cubic symmetry, a J=5/2 multiplet is split into 
a doublet of I\ symmetry and a quartet of Tg symmetry [13]. From the inelastic neutron 
scattering data of reference 9, we know that the doublet has the lower energy in CeCu 2 Si 2 
with the excited quartet about 36 meV higher (i. e. A C ef=3Q meV 360 K). These states 
can be expressed as a linear combination of the J=5/2 states as shown in reference 13, and 
our hybridization matrix elements can be written as 
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V Ta a(k) = J2 c r am V ma (k), (la) 

m=-§ 

where a labels the degenerate states for a given representation T and the expansion coef- 
ficients cram are given in reference 13. The mixing matrix element V mcr has the Coqblin- 
Schrieffer form 

V ma (k) = (-ifaV k] l—^—Y* m _„(k), (lb) 
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where cr=±l is the spin index. 

We build upon the formalism of Read and Newns[14] and calculate the partition function 
as a functional integral over Grassmann and complex fields. The appropriate action can 
ultimately be constructed from the Hamiltonian: 



ka kTa 

+ J2 ^c{k)c\j UqTc s^ + c.c. 

kqTaa 

E tX <l [E fl + qrJkVa + { N m+qr :S kY ~ N ^~^ 



(2) 

Ta " T v 7 J 

where C£ CT and fg ra destroy, respectively a conduction electron of crystal momentum k and 
spin a and an electron in a multiplet state with crystal field quantum numbers T and a. £g is 
the quadratic conduction electron dispersion. Ep represents the bare crystal field energies, 
N r is the degeneracy of the T multiplet, and qr=l/N r . The last term is the constraint 
of unit occupancy at each Ce site, and all energies are measured relative to the chemical 
potential for the free conduction electrons. 

Following Read and Newns, Sg? is the amplitude of a complex field which, at mean- 
field level, renormalizes the hybridization strength, Vr a - One can relate sr to Coleman's 
slave Boson via b r =sre ( - te \ where the phase, 9, can be absorbed into the Fermion via a 
gauge transformation. Notice the rescaled quantities Vr aCT =v^rVrw and % r =s^ r / \/Nr, 
where Vr a and sr are assumed of order one. As elucidated in reference 14, by choosing this 
so-called radial gauge, the theory is not plagued by unphysical infrared divergences. 

At the mean-field level, we assume the s G and A Q fields are uniform in space and time, 
integrate out the Fermionic degrees of freedom and require the resulting free energy be 
stationary with respect to s Q and A G . For simplicity, we have assumed that s r=5 /v / iVr, 
where s Q is the same for both crystal field multiplets. We solve the resulting equations, along 
with an equation for the quasiparticle chemical potential //, self-consistently. See Table 1 
for a set of mean- field parameters that should be applicable to CeCu 2 Si 2 . 
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The quasiparticle states can of course be written as a combination of plane waves and 
crystal field states: 



|Vw(£) >= A n (k) \ka>-Y: ^^A Ta > 

Fa e r - E n (k) 



(6a) 



where er=Er+iA Q are the energies of the shifted crystal field levels and E n (k) is a quasipar- 
ticle energy for the nth band. The normalization factor 



A n {kf 



1 + E 



^ (e r - E n {k)f 



-i 



(66) 



is highly anisotropic. 

It is important to note that the work of Zhang and T. K. Lee, even in the presence of 
anisotropic hybridization, is based on overall spherical symmetry at the cerium ions. Thus, 
their normalization function is isotropic. In our case, the constraints of cubic symmetry 
lead to the strong anisotropies. Specifically, the hybridization with localized states of T 7 

— * 

symmetry vanishes along special directions in k space, resulting in a sharp peak in Ai(k) 
which would actually diverge if the hybridization with r§ states were not nonzero. This 
enhanced normalization affects quasiparticle interactions significantly, as will be discussed 
below. 

Due to the two Bosonic fields, s Q and A c , the dressed Bosonic propagator is conveniently 
calculated from a matrix Dyson's equation D^r^ 1 (q)=D ( -°\ br ^ 1 -Il abir (q). The bare prop- 
agator, 

D (o) 

a b t r, is of order 1/N r , and taking only closed Fermion loops for the self-energy 
insures that the dressed propagator is also of order 1/N r . Allowing the Bose fields to fluc- 
tuate away from the mean-field values and expanding the free energy up to terms quadratic 
in these fluctuations, we can calculate the self-energy U ab ^(q), which in the static limit has 
the following general susceptibility-like structure: 

— * — * 

The effective "matrix element", M ab) r(k,k'), is a complicated anisotropic function, and Q is 
a reciprocal lattice vector. 



The principal value nature of this Brillouin zone sum makes it challenging. We proceed 
by borrowing a numerical technique from electronic structure theory for calculating suscep- 
tibilities: the analytic tetrahedron method [15]. In this procedure we break up the Brillouin 
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zone into a large number of tetrahedra, over which the integral of l/(Ei(fc')-Ei(A;)) can be 
performed essentially analytically. We assume the matrix element is constant inside a given 
tetrahedron, and since M a b t r(k', k) is a complicated function of k and k', this produces fluc- 
tuations in the result as a function of mesh size. This is the most time consuming part of 
the calculation. Our final results for the effective quasiparticle interactions will be presented 
as an average over many different self-energy meshes and will include a standard deviation 
reflecting the fluctuations involved. We wish to emphasize, however, that this "matrix el- 
ement effect" is a well known problem in electronic structure calculations and the current 
approach represents the state of the art [16]. 

The particle-particle scattering amplitude in the even-parity channel is shown in figure 
1. We believe even-parity pairing in CeCu 2 Si 2 is suggested by the following experimental 
evidence: the necessity of strong Pauli limiting to describe the low temperature behavior 
of the upper critical field H c2 (T) (odd-parity would probably produce only a small Pauli 
limiting) [17]; appearance of a large DC-Josephson current through a CeCu 2 Si 2 /Al weak- 
link. (At the simplest level, Al is a singlet BCS superconductor and it is difficult to get a 
Josephson current between opposite parity superconductors.) [18]; and the strong reduction 
of the 29 Si and 63 Cu Knight shifts below T c [19]. 

Because of computational complexity, we confine our calculation to the Fermi surface. 
This weak-coupling calculation means, of course, that we miss any evidence of a pairing state 
with unconventional frequency dependence (a topic which has grown in interest recently). 
We also miss any strong-coupling renormalization of the effective interaction. This should 
not jeopardize, however, our results for instabilities in the E g and T lfl pairing channels, 
since they depend largely on the underlying symmetries of the model and would be material 
dependent. 

For simplicity, we have assumed a spherical Fermi surface. For very small fillings of 
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the lowest quasiparticle band (n toto 2=n c +ii/^Cl) the Fermi surface is indeed truly spherical. 
Only for a nearly filled band (n tota |« 2) should the surface be significantly warped, reflecting 
the Ti symmetry of the problem. We take a filling (n tota i—1.5) that gives a Fermi surface 
only mildly distorted from a sphere and yet still leads to a large quasiparticle density of 
states. 

We have the following form for the particle-particle scattering amplitude (see reference 

8): 

r (k,k') = 1 - £ ^t(^r,,(^) 

xV r *Wa'(k')Vraa'(k') x A 1 (k F ) 2 A 1 (k' F ) 2 x 

~ D {o) jD (°)( :s ° r i S ° T '\]\ 1 l 1 1 (7) 

I s sX \T or + T or J\ Li + F(k> - k) + 1 + F(k> + k)\ ' 1 j 

where T or is an effective Kondo scale for the T crystal field multiplet (T or =er-/u), D s (°)=- 
s r 7 2 /2T or7 , D sA ^=is or7 , and D A ^=0. Recall that A 1 (k) is the anisotropic quasiparticle 
normalization. The function F(A? ± k) has a complicated dependence on the Bosonic self- 
energy. 

For our weak coupling calculation it remains to average the scattering amplitude over 
the Fermi surface 

< r > a = J £ / ^€(k')To(k\ k)Uh (8) 

where the (f> a (k) are the so-called cubic harmonics, i. e. functions that transform under the 
representations of the group O^. The subscript a labels these representations (Ai 9 ,A 29 ,E g , 
Ti g ,T 2fl ). A superconducting instability in the "a-th" pairing channel is signaled by a neg- 
ative value of the corresponding average. 

— * — * 

We now present our main results. First, it is instructive to consider the limit F(k'±k)—> 
(see Table 1). This leaves us with only the local interactions between the quasiparticles, 
which have a large repulsive strength in the A lfl (s-wave like) channel and are enhanced 
by the strong anisotropy of the normalization function. Furthermore, a combination of 
the normalization and the hybridization matrix element causes the local interactions to be 
attractive in the E^d^^) channel and repulsive but small in the T\ g and T 2ff channels. 
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Our results contrast with those of Zhang and Lee[8] . They find repulsive local interactions 
in the s,d, and g wave pairing channels that are all approximately the same size. We 
emphasize that our new results stem from the lowering of the symmetry from spherical to 
cubic. 

The last step is to include the full bosonic self-energy through the function F(k'±k), 
which is strongly anisotropic due to its dependence on the effective matrix element, M^r- 
This anisotropy gives rise to an averaged scattering amplitude (equation 8) that varies with 
the self-energy mesh size. We have thus averaged the results of equation 8 over eleven 
different self-energy meshes, and the results for the effective quasiparticle interactions are: 

\ Alg = ~ < r o > Alg /T r 7 = -1.70 ± 0.04 

X Eg = ~ < r o > Eg /T orr = 0.133 ± 0.10 

X Tlg = -<T > Tlg /T oT7 = 0.231 ± 0.18 

= - < r o > T2g /T oT7 = -0.616 ± 0.45. 

The fluctuations may be large because of the matrix elements in the Bosonic self-energy, 
but a superconducting instability in the E g and T l3 channels is still clear. 

We wish to remark that a strong-coupling calculation, including the dynamics of the 
exchanged boson, is currently intractable with this procedure. However, we don't expect it 
to modify the qualitative aspects of our stability analysis since the quasiparticles interactions 
at T c are renormalized by factors of 1/(1 + A) where A arises from the frequency dependence 
of the quasiparticle self energy. Since the quasiparticle self energy in heavy fermion materials 
is expected to be largely momentum independent, this shouldn't affect the ordering of the 
Ar couplings. Furthermore, because of the complexity of reliably computing the self-energy, 
pushing even the static calculation to order 1/N 2 also appears unrealistic. The anisotropies 
due to the cubic symmetry have increased the complexity of the problem more than one 
might have expected. We do imagine, however, studying other Ce based systems such as 
CeAl 3 at the 1/N level [20]. 
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Figure: Quasiparticle scattering amplitude in the pseudospin singlet channel. The full 
lines represent quasiparticles at the Fermi surface; the dashed line is the dressed boson 
propagator. The large circles are bare, anisotropic vertices that connect quasiparticles with 
the bosons. 
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TABLES 

TABLE I. Fermi Surface Average of the Local "Hard-Core" Quasiparticle Scattering Ampli- 
tude, "a" labels the irreducible representations of the group Cv "local" means the bosonic 
self-energy was not included. The results from Zhang and Lee[8] are also presented for comparison. 
It is clear that cubic symmetry yields different quasiparticle interactions than spherical symmetry. 
The mean- field parameters used in our calculations are: Kondo temperature (T O 7=1.30 meV); 
crystal field splitting (Acef=0-036 eV); hybridization renormalization (s o =0. 1424094); and bare 
hybridization strength (V o =0.65 eV). 
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Cubic Crystal Fields 


Spherical Symmetry (from [8]) 


a 


(r o ) a /T r7 


1 


(r )i/T 


Ai 3 


3.79 


1=0 (s- wave) 


1/6 


Eg 


-0.019 


1=2 (d- wave) 


4/21 


Ti 9 


0.130 


1=4 (g- wave) 


1/7 


T 25 


0.133 
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